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ABSTRACT 

The H2O Southern Galactic Plane Survey (HOPS) has mapped a 100 degree strip of 
the Galactic plane (-70° > I > 30°, \b\ < 0.5°) using the 22-m Mopra antenna at 
12-mm wavelengths. Observations were conducted in on-the-fly mode using the Mo- 
pra spectrometer (MOPS), targeting water masers, thermal molecular emission and 
radio-recombination lines. Foremost among the thermal lines are the 23 GHz transi- 
tions of NH 3 J,K = (1,1) and (2,2), which trace the densest parts of molecular clouds 
(n > 10 4 cm~ 3 ). In this paper we present the NH 3 (1,1) and (2,2) data, which have 
rS , a resolution of 2 arcmin and cover a velocity range of ±200kms _1 . The median sen- 

■ sitivity of the NH3 data-cubes is <7T mb = 0.20 ± 0.06 K. For the (1,1) transition this 

sensitivity equates to a 3.2 kpc distance limit for detecting a 20 K, 400 M Q cloud at 
the 5cr level. Similar clouds of mass 5,000 M Q would be detected as far as the Galactic 
centre, while 30,000 M Q clouds would be seen across the Galaxy. We have developed an 
automatic emission finding procedure based on the ATNF duchamp software and have 
used it to create a new catalogue of 669 dense molecular clouds. The catalogue is 100 
percent complete at the 5a detection limit (T m b = 1.0 K). A preliminary analysis of 
the ensemble cloud properties suggest that the near kinematic distances are favoured. 
The cloud positions are consistent with current models of the Galaxy containing a 
long bar. Combined with other Galactic plane surveys this new molecular-line dataset 
constitutes a key tool for examining Galactic structure and evolution. Data-cubes, 
spectra and catalogues are available to the community via the HOPS website. 

Key words: starsdormation, ISM:evolution, radio lines:ISM, Galaxy: structure, sur- 
veys, stars:early type 



1 INTRODUCTION 
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HOPS (H2O Southern Galactic Plane Survey) is a project 
utilising the Mopra radio telescop^B to simultaneously map 
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spectral-line emission along the southern Galactic plane 
across the full 12- mm band (frequ encies of 19.5 to 27 .5 GHz). 
Since the survey began in 2007 ijWalsh et al.ll2008l ) HOPS 
has mapped 100 square degrees of the Galactic plane, from 
/ = 290°, continuing through the Galactic centre to I = 30° 
and with Galactic latitude |6| 0.5°. The aim of the sur- 
vey is to provide an untargeted census of 22.235 GHz H2O 
(616— 523) masers and thermal line emission towards the in- 
ner Galaxy. Observations were completed in 2010 and the 
survey properties, observing parameters , data reduct i on and 
H2O maser catalogue are described in IWalsh et al.l (|201ll ) 
(hereafter Paper I). 

The primary thermal lines targeted by HOPS are those 
of ammonia (NH3), whose utility as a molecular thermome- 
ter in a broad range of environments is unsu rpassed. With 
an ef fective critical density of ~ 10 cm" " 3 |Ho fc Townesl 
1983), NH3 traces dense molecular gas and it is often asso- 
ciated with the hot molecular core phase of high-mass sta r 
formation (e.g.. lLongmore et al. 1 12001 iMorgan et aT1l2010h . 
where it exhibits temperatures in excess of 30 K. NH3 is ex- 
cited in gas with k inetic temperatures greater than ~ 5K 
l|Pickett et al.| [l998) and is also found associated with cool 
(T < 10 K) dense clouds. Such regions are too cold for 
more common gas tracers, like CO, to remain in the gas 
phase. Instead they are fro zen out onto the surfaces of dust 
grains (jBergin et al.ll2006l ). The J,K=(l,l) inversion tran- 
sition exhibits prominent hyperfine structure, which can be 
used to infer the optical depth of the transition. In clouds 
forming high- mass stars (^ 8Mq) and under optically thin 
conditions, the peak brightness of the four groups of satel- 
lite lines is approxim ately half that of the central group 
|Rvdbeck et al.|[l977T ). Comparison of the (1,1) and higher 
J,K inversion transitions can be used to estimate the rota- 
tional temperature of the gas. 

In this second paper we present the HOPS NH3 (1,1) 
and NH3 (2,2) datasets and the automatic finding procedure 
used to create catalogues of emission. We illustrate the basic 
properties of the catalogues, which form the basis for fur- 
ther analysis. A third paper (Longmore et al. in prep.) will 
published the properties of the catalogue derived by fitting 
the NH3 (1,1) and NH3 (2,2) spectra with model line profiles 
(i.e. temperature, density, mass and evolutionary state). 



2 OBSERVATIONS AND DATA REDUCTION 

The HOPS observa tions and survey design are described in 
IWalsh et alj (|2008T l and in Paper I. For convenience a brief 
summary is provided here. 

2.1 Observations 

Observations were conducted using the 22-m Mopra radio- 
telescope situated at latitude 31:16:04 south, longitude 
149:05:58 east and at an elevation of 850 meters above sea 
level. Data were recorded over four summer seasons during 
the years 2007-2010 using a single-element 12-mm receiver. 

The digitised signal from the receiver was fed into 
the Mopra Spectrometer (MOPS) backend, which has an 
8.3 GHz total bandwidth split into four overlapping sub- 
bands, each 2.2 GHz wide. During the HOPS observa- 
tions the spectrometer was configured in 'zoom' mode, 



whereby each of the sub-bands contained four 137.5 MHz 
zoom bands. Up to sixteen zoom bands may be recorded 
to disk, each split into 4096 channels. A single zoom 
band covered both the J,K=(1,1) and (2,2) NH 3 in- 
version transitions, which have line-centre rest frequen- 
cies of 23.6944803 GHz and 23.7226336 GHz, respectively 
|Pickett et al.| [l9981. This frequency setup resulted in a ve- 
locity resolution of ~ 0.42 km s -1 per channel, spanning a 
velocity range of 1739 km s -1 . 

At 23.7 GHz the full-with half- maximum (FWHM) o f 
the Mopra beam is 2.0 arcmin (|Urauhart et all bOlOp . 
meaning that the 100 square degree target area would re- 
quire 90,000 pointings to create a fully sampled point-map. 
Based on experien ce gained during the pilot observations 
l|Walsh et al.ll2008l 1 the survey area was divided up into 400 
0.5° x 0.5° maps in Galactic coordinates. The telescope was 
driven in on-the-fly (OTF) mode, raster-scanning in either 
Galactic / or b and recording spectra every six seconds. At 
the end of each row spectra were taken of an emission-free 
reference position situated away from the Galactic plane. 
Scan rows were offset by half a beam FWHM and a scan 
rate was chosen to ensure Nyquist sampling at the high- 
est frequency. Each map was observed twice by scanning 
in orthogonal I and b directions the results were averaged 
together to minimise observing artifacts. 



2.2 Data reduction 

Raw data from the telescope were written into rpfits for- 
mat files consisting of lists of spectra with associated time- 
stamps and coordinate information. The ATNF lived att^] 
software package was used to apply a bandpass correction 
by forming a quotient between the reference spectrum and 
individual spectra in its associated scan row. The software 
was configured to fit and subtract a first order polynomial 
from the bandpass, excluding any strong lines and 150 noisy 
channels at each edge. 

The gridzilla package was used to resample the band- 
pass corrected spectra onto a regular coordinate grid. The 
software also performed interpolation in velocity space to 
convert measured topocentric frequency channels into LSR 
velocity. At this stage spectra with system temperatures 
T sys > 120 K were discarded. Software settings were cho- 
sen to produce data cubes with pixel scales of 30 x 30 arcsec 
and a resolution of 2 arcmin. The final data products were 
1 x 10 deg 2 cubes with 3896 usable velocity channels. For 
each spectral line within this pass-band the data has been 
cropped to ±200 km s" 1 centred on the line centre veloc- 
ity (determined from the line rest frequency), which corre- 
sponds to the normal velocity ra nge for molecular material 
in the Galaxy (|Dame et al.ll200lT ). 

A small number of spectra (~ 10 percent) were affected 
by a baseline ripple or a depressed zero level. We fitted ad- 
ditional polynomial baselines of order 3-5 to the line- free 
channels of these data. See Appendix 1X1 for more detail. 



2 http://www.atnf.csiro.au/people/mcalabre/livedata.html 
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Figure 1. Map of the peak NH3 (1,1) signal-to-noise (S/N) across the HOPS survey area. The map was made by smoothing the 
S/N cube spatially (final beam FWHM = 2.5 arcmin) and in velocity (hanning- window = 5) and measuring the peak brightness in 
each spectrum. An approximate brightness temperature scale may be found by multiplying by the average root-mean-squared noise 
temperature (ot b ) = 0-20 K. 
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Figure 2. Map of the peak NH3 (2,2) signal-to-noise (S/N) across the HOPS survey area. 
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Figure 3. The top panel shows the sum of the NH3 (1,1) inte- 
grated intensity over the survey region as a function of Galac- 
tic longitude in four-degree bins. The longitude bin covering 
the Central Molecular Zone (CMZ, -2° < I < 2°) peaks at 
J T mb dv = 4.65 X 10 5 Kkms - 1 1 and the y-axis scale has been 
truncated so as to show the remainder of the emission. The inset 
shows how the plot appears scaled to the full J T mb dv range. 
In the bottom panel is a similar plot of J T mb dv versus Galactic 
latitude. The filled grey histogram includes all emission, while the 
hatched histogram includes only emission outside the CMZ. 



2.3 Temperature scale and uncertainty 

Data from Mopra are calibrated in brightness temperature 
units (K) on the T£ scale, i.e., they are corrected for radia- 
tive loss and rearward scattering (although not atmospheric 



atten uation at 12-mm wavel engths - see iKutner fc Ulichl 
Il98ll for a review and also lLadd et al.l 120051 ). It is de- 
sirable to convert this onto a telescope-independent main 
beam brightness temperature scale (T mb ), which is also cor- 
rected for forward scattering. The measured T mb of a source 
would equal the true brigh tness temperatu r e if th e source 
just filled the main beam. lUrquhart et al.l (|2010t) charac- 
terised the telescope efficiencies between 17 GHz and 49 GHz 
and we divided the HOPS data by the main-beam effi- 
ciency, ») m b, given therein. For both NH3 (1,1) and (2,2) data 
?7mb = 0.571. The uncertainty on the T mb scale for HOPS is 
~ 30 percent (|Walsh et al.ll201ll ). 
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Figure 4. Plots of the brightness ratio T mb n ,i)/T, 



mb, (2,2) 



for 



every emitting voxel (l-b-v data element) with T mb > 3<r in both 
HOPS NH3 datasets. The hatched histogram includes only voxels 
from the CMZ (|/| < 2°) while the solid histogram includes the 
remaining emission from the outer Galaxy. 



3 THE NH 3 DATA 

In contrast to other large molecular line surveys of the 
Galactic pla ne, which use rela t ively abundant CO isotopo- 
logues (e.g. iDame et all l200ll , Ijackson et ail l200rj ) HOPS 
has focused on the inversion-rotation transitions of NH3. 
Detections derive from the densest parts of giant molecular 
clouds where gas is condensed into cool clumps, or in the 
case of warm gas (T > 3 K) surrounding ce ntrally heated 
star-forming regions (e.g.. iPurcell et aT1l2009t ). 



3.1 NH3 emission properties 

Figures [T] and [2] present peak signal-to-noise (S/N) maps 
covering the NH3 (1,1) and (2,2) passbands, respectively 
(I^lsrI < 200 kms -1 ). Such maps are the easiest way to 
visualise the distribution of emission on the sky without be- 
coming confused by noise-related artifacts. The maps were 
constructed by smoothing a S/N cube to a spatial resolution 
of 2.5 arcmin and then hanning-smoothing in the spectral di- 
mension using a filter width of five channels (~ 2.0 kms -1 ). 
For any given spatial pixel the brightest spectral channel was 
recorded on the map. The S/N cube was generated during 
the emission-finding process, which is explained in Section^ 



3.1.1 Galactic distribution of emission 

The NH3 (1,1) map (Figure [TJ reveals a multitude of dis- 
tinct cloudfl between -30° < I < 30°. The density 
of clouds falls off rapidly in the southern Galaxy be- 
yond I < —30°. Several well known star-formation com- 
plexes are detected as bright peaks including G33 3.0-0.4 



plexes are detecte d as brignt p 
jBains et al.ll2006l, I Wong et all [2 



120081. lLo et al. I2009D . W43 



(G29.9.lDame et al.lll98fj,lNguven Luong et al.ll201ll) . G305 
( Clark fc Porteill2004l , [Hindson et alj|2010l) and the 'Nessie ' 
filament at I = -21.5° , b = -0.48° i| Jackson et alj|2010h . 



3 We refer colloquially to any contiguous region of emission as 
a 'cloud' without implying membership of any object category 
based on size-scale. 
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Figure 5. Detail of the longitude-velocity diagram showing 
the Galactic centre reg ion (top). NH3 (1, 1) emission is shown 
in colourscale and the iDame et all ll200ll ) CO (1-0) contours 
are overplotted. While the overall emission morphology of both 
molecules is similar, the brightest peaks are offset in / and v. 
The emission seen at v < — 170 km s is due to broad linewidth 
NH3 (2,2) emission encroaching on the NH3 (1,1) pass-band. The 
bottom panel shows the NH3 (1,1) integrated intensity map for 
the same I range. 



Emission from the Galactic centre region is most promi- 
nent between —2° < I < 2° and although symmetric 
in overall morphology, is significantly brighter at positive 
longitudes. Gas in this region, known in the literature as 
the Central Molecular Zone (CMZ), is characterised by 
high temperatures and densities, and by large velocity dis- 
persions as it flows into the central 200 pc of the Galaxy 
ijMorris fc Serabvnlll996h . 

Figure [3] illustrates the distribution of NH3 (1,1) emis- 
sion as histogram. The top panel shows the total integrated 
intensity (summed over all v and b) as a function of Galac- 
tic longitude in 4° bins. The CMZ contains 80.6 percent 
of the detected NH3 (1,1) emission and dominates the plot, 
which is scaled to show the less intense emission in the outer 
Galaxy. The bottom panel shows the same plot as a func- 
tion of Galactic latitude in 3 arcmin bins. The filled grey 
histogram incorporates all emission in the HOPS target re- 
gion. The broad and intense peak at b — —0.075° is solely 
due to emission from the CMZ. However, with the Galactic 
centre removed (hatched histogram), the distribution with 
b is relatively flat. 

It is interesting to note tha t the Galactic dist ribution 
of H 2 masers found in HOPS (| Walsh et alJl201ll ) is very 
different to that of the dense molecular gas traced by NH3. 
The masers show no evidence of a longitude peak at the 
Galactic centre, but are clustered around the mid-plane of 
the Galactic disk with an angular scale-height of (9fwhm = 
0.60°. The implications for star-forming environments in the 
CMZ compared to the spiral arms will be discussed in a 
future paper (Longmore et al., in prep). 
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Figure 6. Integrated NH3 (1,1) spectrum for the Galactic centre 
(top panel), the region around the G3. 34 + 0. 42 cloud, also known 
as Bania's Clump 2 (middle panel), and the outer Galaxy (bot- 
tom panel). Spectra were integrated over the full latitude range 
(±0.5°). 



3.1.2 NH 3 (1,1) versus NH 3 (2,2) emission 

Only one-third of the NH3 (1,1) emission regions detected 
above 3a have counterparts in the NH3 (2,2) data, de- 
spite the data having similar sensitivities. Most clouds with 
T mbi n t i) >1 K are detected in the (2,2) data, although there 
are exceptions (for example the cloud at I = — 2.5°, b — 
0.33°). Where emission is detected in both lines, the bright- 



ness ratio T, 



mb, (l,l)/r m b, (2,2) 

over the HOPS area. Clouds with low ratios or no detected 
NH3 (2,2) emission tend to be associated with extinction 
features in the GLIMPSE0 8 urn infrared images and are 
excellent candidates for the cold and dense precursors to 
clusters of high-mass stars. Strikingly, the morphology of the 
CMZ is almost identical in both datasets and the brightness- 
temperature ratio is uniformly close to one. Figure 2] illus- 
trates this difference graphically by plotting the distribu- 
tion of brightness-temperature ratios between every com- 
mon emitting voxel (3-dimensional data elements in I, b & 
v) with T m b > 3cr. The tall, narrow histogram (hatched) 
contains voxels drawn from the CMZ only and has a median 
of 1.2 and a width of a = 0.19. By comparison the data 
from the outer Galaxy has a median of 1.6 and a = 0.50. 



varies between ~ 0.2 and ~ 4.2 



4 Galactic Legacy Infrared Mid-Plane Survey Extraordinaire 
http://www.astro.wisc.edu/sirtf/ 
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Lower ratios in the CMZ are consistent with warmer exci- 
tation conditions, but may also be due in part to the effects 
of high optical depths or line saturation. 

A small proportion of NH3 (2,2) emission is not asso- 
ciated with NH3 (1,1) emission at the same location. The 
J,K = (1,1) transition is more easily excited than the (2,2), 
hence such NH3 (2,2) detections are most likely artifacts. 
We analyse the completeness and spurious source counts in 
Sections [5] and 15.31 



3.1.3 Position-velocity Diagram 

Spectral lines originating from the Galactic centre region 
have extremely broad profiles (> 50kms _1 ) and merge to- 
gether into a single bright region of emission with a com- 
plex morphology. Figure \E\ (top) is an l-v diagram show- 
ing the inner eight degrees of the Galactic plane and the 
prominent CMZ. To make the map, voxels without detected 
emission were set to zero and the data were summed in 
Galactic latitude. The corresponding integrated intensity 
map is presented in the bottom panel. A significant ve- 
locity gradient exists across the emission , which has also 
been observed in the Dame et al. J200 J) CO (1-0) data. 
iRodriguez-Fernandez fc Combe si (|2008l ) have modelled the 
CMZ as the gas response to a short 'nuclear bar', which it- 
self is embedded within a longer Galactic bar, thought to be 
inclined towards us at an angle of 25° — 35° , with the nearer 
end at positive Galactic longitudes. The large broad-line re- 
ion of emission at I — 3.1° is known as Bania's Clump 2 



Banialfl977l ) and has been theorised to lie at the closest 
end of innermost elonga ted x-1 orbit around the Galactic 
centre |Ballv et alj|2010h. The CO (1-0 ) contours from the 
12 arcmin resolution iDame et alj (|2001) survey are overlaid 
on the l-v diagram for comparison. The overall morphology 
of the emission is similar, however, the 2 arcmin resolution 
NH3 data shows more structure, both spatially and spec- 
trally. Galactic CO (1-0) emission is optically thick com- 
pared to NH3, which probes the kinematics of the molecu- 
lar clouds at greater depths. NH3 detections in HOPS are 
severely distance limited (see Section I5.3.3|l and detecting 
emission from beyond the Galactic centre is difficult, even 
in the highly excited central Galaxy. 

3.1.4 Integrated spectrum 

The NH3 (1,1) spectrum is characterised by five groups of 
emission lines seperated ~ lOkms" 1 . Under optically thin 
conditions the central group is approximately twice the 
brightness of the satellites. Figure[6]presents NH3 (1,1) spec- 
tra summed over all b and selected I ranges. In the top 
panel is shown the integrated spectrum of the CMZ be- 
tween —2° < / < 2°. Individual spectra (and line-groups 
within spectra) are blended into a single emission feature 
spanning -leOkms" 1 < ulsr < leOkms" 1 . The middle 
panel shows the integrated spectrum of the region domi- 
nated by Bania's clump 2 (G3.34 + 0.42 in this work), which 
has the broadest linewidths outside of the CMZ. Emission at 
VhSR < ~ 170 kms -1 is due to broad NH3 (2,2) lines imping- 
ing on the NH3 (1,1) bandpass. The spectrum in the bottom 
panel is constructed from the outer Galaxy data only and 
contains multiple distinct velocity features deriving from in- 
dividual regions of emission. The peaks at approximately 



— SOkms^ 1 and +40kms _1 are likely associated with the 
Galactic ring and Scutum-Centarus spiral arm, respectively 
(see Section Some of the complexity evident in the en- 
semble spectrum may also be attributed to the summation 
of satellite line-groups. 



3.2 Noise characteristics 

Figure [7] shows an image of the root-mean-squared (RMS) 
noise temperature or mb over the whole survey region. The 
image was constructed by measuring the RMS noise along 
the spectral dimension of the NH3 (1,1) cube after mask- 
ing off all channels containing line-emission. We are con- 
fident that the noise map does not contain any contribu- 
tion from real emission as contamination from the broad- 
linewidth CMZ is absent, and this represents the worst-case 
scenario. 

The defining feature of Figure [7] is the chequered pat- 
tern corresponding to individual 0.5 x 0.5 degree maps. The 
variable noise level between maps is due to temporal changes 
in the conditions encountered while observing. Within each 
map, rapidly fluctuating cloud cover and flagged or missing 
spectra lead to high- noise 'hot-spots'. More gradual changes 
in observing conditions (e.g., humidity and elevation) mani- 
fest themselves as significant noise variations between scan- 
rows in I or b. Automatically finding line-emission in such 
an inhomogeneous data-cube poses a distinct challenge. 

Histogram plots of the noise temperature CT mb for all 
spatial pixels in the HOPS NH3 (1,1) and (2,2) data are 
plotted in Figure [8] Both datasets have Gaussian shaped 
distributions, except for a small number of pixels which oc- 
cupy a tail extending to higher noise temperatures (0.71 % 
above 3a). The median noise temperature of the J,K = (1,1) 
data is higher than the (2,2): <JT mb (1 y = 0.201 K compared 
to <TT mb (2 9) = 0.186 K, although the Gaussian part of both 
distributions have similar widths at FWHMj^!) = 0.056 K 
and FWHM(2.2) = 0.053 K. Upon inspection of the data we 
find that the high-noise pixels in the tail arise from ~ 10° 
maps observed during poor weather and from the edges of 
the maps where fewer raw spectra contribute to individual 
pixels. The difference between the NH3 (1,1) and (2,2) distri- 
butions stems from a single low-noise NH3 (2,2) map which 
had additional data added into the processing pipeline. No 
emission was detected in this map, which lies in the far 
southern Galactic plane at I = —60.25, b — 0.25. 



4 SOURCE FINDING AND MEASUREMENT 

The full HOPS dataset covers 100 square degrees of the 
sky, equivalent to over 90,000 overlapping Gaussian beams 
at 23.7 GHz. For each position we obtain a spectrum span- 
ning a 400kms _1 velocity range, or 931 channels. Manually 
identifying emission in such a large volume of data would 
be prohibitively tedious and likely error-prone, so an au- 
tomatic emission finding procedure was written for HOPS. 
At the core of the method is the ATNF's DUCHAMF0 soft- 
ware (Whiting 2012), which has been developed to automat- 
ically detect emission in three-dimensional data, duchamp 



http: / /www. atnf.csiro.au/people/Matthew.Whiting/Duchamp/ 
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Figure 7. Image of the root-mean-squared noise temperature in T mb units for the entire HOPS region made from the NH3 (1,1) data. 
Individual observing blocks are visible as 0.5° chequered squares, while scanning artifacts present as lines of high or low noise in Galactic 
£ and b. The median root-mean-squared noise is 0.20 K± 0.05 K, where the uncertainty is the FWHM of the total noise distribution (also 
see Figure \E\ 
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Figure 8. Histograms of the RMS noise temperature for all (I, b) 
spatial pixels in the NH3 (1,1) and (2,2) data (grey-shaded and 
hatched histograms, respectively). The noise temperatures were 
determined from the line-free channels in the spectra. Inset is the 
same histogram, but with a logarithmic scaling on the Y-axis to 
enhance the view of the high-noise tail. 
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Figure 9. Example signal-to-noise spectrum (grey line) extracted 
from a single pixel with the wavelet-reconstructed spectrum over- 
laid (red-line). DuCHAMP does an excellent job of suppressing the 
noise allowing the use of a low threshold: 0.8cr in this case. 

is specifically designed for the case of a small number of sig- 
nal voxels within a large amount of noise and is perfectly 
suited to the wide bandwidth data-cubes produced by the 
Mopra spectrometer. 

4.1 Overview of the source finding procedure 

DUCHAMP searches a cube by applying a single threshold, 
either a flux value or a signal-to-noise ratio, for the whole 
dataset, and so is thus best applied to data with uniform 
noise. To compensate for the variable noise in the NH3 data, 
the HOPS emission finding procedure implements a two-pass 
solution. Firstly DUCHAMP is run on a smoothed version of 
the entire cube down to a global 3a cutoff level. A mask 
data-cube is produced which is used to blank the emission 



in the original data. Some of this 'emission' identified in the 
first pass may in fact be regions of high RMS noise and 
conversely some real low-level emission may have escaped 
detection. However, this likely does not matter as the goal 
of this step is to create a cube containing mostly line-free 
noise. For every spatial pixel in the blanked cube the stan- 
dard deviation (o"T mb ) is estimated along the spectral axis 
to create a map of background noise temperature like that 
shown in Figure [7] In practice we calculate the median ab- 
solute deviation from the median (MADFM), which for a 
dataset X = xi, X2 . . . Xj. . . x n is given by: 

a — K. median (|xi — median(X)|), (1) 

i.e., the median of the deviations from the median value. 
In Equation [T] K is a constant scale factor which depends 
on the distribution. For normally distributed data K « 1.48. 
The MADFM statistic is largely robust to the presence of a 
few channels much brighter than the noise. This is especially 
important when estimating the noise in spectra containing 
very broad lines, such as those observed towards the Galactic 
centre. 

Each spectral plane of the original input cube is then di- 
vided by the noise map to make a signal-to-noise cube with 
homogeneous noise properties. DUCHAMP offers the option of 
reconstructing data-cubes using the a trous wavelet method 
prior to running the finder. A thorough description of the 
procedure may be found in IStarck fc Murtaghl |l994). The 
reconstruction is very effective at suppressing noise in the 
cube, allowing the user to search reliably to fainter levels 
and reducing the number of spurious detections. We chose 
to reconstruct the NH3 data in 3D-mode, meaning that the 
wavelet filters can distinguish between narrow noise spikes 
confined to single channels and small regions of emission 
spanning l-b-v space. Figure [9] shows an example NH3 (1,1) 
spectrum with the reconstructed version overplotted. A sec- 
ond pass of DUCHAMP is run on the reconstructed cube to 
construct the final list of emission sources. A detailed de- 
scription of the DUCHAMP inputs is presented in Appendix iBl 

4.2 Cloud measurements 

The properties of NH3 clouds were measured directly from 
the NH3 data cubes using the 3D voxel-based emission 
masks produced by the source finder. Two dimensional pixel- 
masks were made by collapsing the 3D masks along the ve- 
locity axis. 

4-2.1 Position, velocity and angular size 

Three position and velocity measurements were made on 
each cloud: centroid, weighted and peak. The centroid po- 
sition (l c , b c , v c ) was measured directly from the 3D voxel 
mask and corresponds to the geometric centre of the cloud 
in l-b-v space. A brightness weighted position (Z w , 6 W , «w) 
was measured by weighting each voxel coordinate with its 
brightness temperature value T m b according to: 

n J n 

= ^ ^ ni T^nb,! / ^ ^ ?mb,i , (2) 
i = l / i = l 

where n is the coordinate axis to be measured. The peak 
position (l p , bp, v p ) is the coordinate of the brightest voxel 
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in the 3D cloud. Noise spikes can introduce significant errors 
in the weighted and peak measurements, especially for weak 
detections, so we performed the coordinate measurements 
on cubes spatially smoothed to a resolution of 2.5 arcmin 
and spectrally smoothed using a hanning window five chan- 
nels wide. Sources with simple morphologies have similar 
positions and velocities in all three measurements compared 
to extended sources in which the measurements differ. 

4-2.2 Angular size and area 

The angular size of the clouds was quantified in two ways. 
Firstly, the solid angle Q subtended by the cloud was mea- 
sured from the 2D pixel mask. The equivalent angular ra- 
dius r c is then defined as the radius of a circular source 
which would subtend an equivalent solid angle on the sky. 
Secondly, the brightness-weighted radius r w was calculated 
with respect to the (l w , 6 W ) position via: 

n / n 

= ^ ^ T% Tmbji / ^ ^ ^mb , i ; (3) 
i=l / i=l 

where r; is the distance between the brightness-weighted 
centre and the i* pixel. For a high signal-to-noise ratio un- 
resolved source this is directly equivalent to the Gaussian 
FWHM/2. 

4-2.3 Integrated intensity and brightness temperature 

The integrated intensity J T m b dv of a cloud is the summed 
brightness of the emitting voxels (in K) times the velocity- 
width Aii of a channel (in kms" 1 ). We have an excellent 
measurement of the average RMS noise over the cloud (see 
Figure [7}, so the uncertainty on J T m b dv is given by: 

5(jT mb dvj = (RMS)AvVN, (4) 

where N is the number of independent spatial pixels sub- 
tended by the cloud. Systematic fluctuations in the spectral 
baselines may also contribute to the error, but these are 
likely insignificant compared to the HOPS sensitivity. 

The peak brightness temperature T m b is simply the 
value of the brightest voxel in the cloud. 

4.3 Example clouds 

Forty-four percent of NH3 (1,1) clouds in HOPS are unre- 
solved (r w < 1 arcmin). Figure [10] presents examples of two 
bright and extended clouds. Both clouds show evidence of 
sub-structure in position and velocity, but are merged into a 
single object by the emission finder. G003. 145+0301+066. 5 
(named for l c +b c +v c ), also known as Bania's Clump 2, 
(Banral ll977l ). is interesting as it is one of the few clouds 
outside of the CMZ to exhibit very broad line profiles (8v > 
50 km s~ 1 ) . The emission is largely confined to two lobes sep- 
arated in b and v, but with components which merge at v ~ 
75 kms" 1 . Physical conditions in G003. 145+0301+066.5 are 
comparable to the CMZ, as the median brightness ratio 
Tmb, (1,1) /Tmb, (2,2) = 1-2 compared to 1.6 for clouds in the 
outer Galaxy. 

G333.183— 0.405— 52 is a well known giant molecular 



cloud containing Hn regions, ho t molecular cores and out- 
flow sources fsee lLo et alj [20091 and references therein). In 
contrast to Bania's Clump 2, the line profiles are well- 
defined. NH3 (1,1) exhibits the classic spectral shape at both 
sampled positions: a central main group of lines with four 
symmetric satellite groups, approximately half as bright as 
the main-group. The spectrum o f NH3 (2,2) also ex hibits 
four groups of satellite lines (see iPickett et alH l998), but 
they are generally too weak to be detected in the HOPS 
data. 

Overlaid on the J Tmbdv maps are the three position 
measurements as described in Section [4]2] The peak voxel is 
denoted by a star, the geometric centre is marked by a '+' 
and the brightness- weighted centre by a ' x ' symbol. Note 
that the peak voxel is not necessarily coincident with the 
peak pixel of the J T^dv map. Geometric and brightness- 
weighted positions are generally in close agreement, even 
in clouds with complex morphology. Indeed, separations 
of more than a beam-width between these two positions 
and the peak voxel position indicates that multiple clumps 
of emission exist within a cloud boundary. The weighted 
method is the most robust, so we adopt the weight values 
from now on we when referring to longitude, latitude and 
fLSR of a source. 



5 THE HOPS NH 3 CATALOGUE 

The emission finder was run on the NH3 (1,1) and (2,2) 
datasets independently to produce the final catalogues of 
clouds. In total 669 NH 3 (1,1) and 248 NH 3 (2,2) clouds were 
found in the HOPS cubes. In the following sections we exam- 
ine the catalogues produced by the source finding routines. 
We aim to quantify the robustness of the duchamp detec- 
tions and investigate the distribution of source properties. 



5.1 NH 3 (1,1) versus (2,2) catalogues 

The initial NH3 (1,1) catalogue contained 687 detections of 
which 18 were flagged as artifacts. Five detections between 
— 2 < I < 2 with «lsr < — WOkms" 1 are in fact bright, 
extended clouds of emission from the NH 3 (2,2) transition. 
The extremely large linewidths close to the Galactic centre 
(50 kms" 1 < Av < 100 kms" 1 ) cause NH3 (2,2) emission to 
spill over into the NH 3 (1,1) bandpass. A further thirteen 
detections are flagged as they arise from artifacts caused by 
the baselining routine. Due to the polynomial used to fit 
the baseline (see Appendix |A")) , the spectra at the edge of 
the bandpass can rise slightly above the noise cutoff for sev- 
eral contiguous velocity channels. If this occurs for enough 
contiguous spatial pixels to pass the detection threshold the 
emission will be classed as a source. However, such sources 
are easy to flag by eye. 

Duchamp originally found 324 individual clumps 
within the NH3 (2,2) data-cube of which 76 are flagged as ar- 
tifacts. Eighteen spurious detections are due to a turned-up 
bandpass-edge, or an intrusion by a broad NH3 (1,1) linew- 
ing. The J,K = (1,1) transition is more easily excited than 
the (2,2), hence NH3 (2,2) clouds are only considered valid 
if associated with (1,1) emission at the same position. Fifty- 
eight low level NH3 (2,2) detections (T m b < 5a, n p i xc i s ^ 12) 
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Figure 10. Examples of two bright clouds in HOPS. For each cloud we show the NH3 (1,1) and (2,2) integrated intensity maps (bottom- 
left and bottom-middle panels, respectively), summed longitude velocity (l-v) diagrams (top-left and top-middle panels) and sample 
spectra from two local peak positions (right panels). The l-v diagrams were made by summing over the Galactic b axis of the cubes. The 
position of the brightest voxel is marked by a star symbol, the geometric centre by a '+' and the brightness-weighted centre by a 'x'. 
See section [4.31 for details. 



were flagged as spurious for this reason. In Section [5.31 be- 
low, we investigate the formal catalogue completeness and 
contamination by spurious sources. 



5.2 Catalogue description 

Tables [2] and [3] present a sample of the catalogue entries 
for the NH3 (1,1) and (2,2) data, respectively. The complete 
catalogue is available in electronic format on the HOPS web- 



sitqj. Columns in the tables are as follows: (1) the cloud 
name constructed from the centroid position of the cloud in 
Galactic longitude, latitude and velocity; (2) the brightness- 
weighted Galactic longitude, (3) latitude and (4) vlsr mea- 
sured from a masked and smoothed version of the original 
cube. Columns (5), (6) and (7) contain the Galactic coordi- 
nates (I & b) and velocity (wlsr) of the brightest voxel in the 
cloud. Columns (8) and (9) tabulate the velocity range over 
which duchamp detects emission within that clump; column 
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Figure 11. Histogram of the number of spurious sources detected 
in an 'empty' NH3 (2,2) data cube (see Section l5,3.1l for details). 
The source finder was run at a range of sigma-limits to charac- 
terise the number of detections as a function of source solid angle. 
Grey bars = 0.7o", red-hatched (left-slanted) bars = 0.8cr and blue- 
hatched (right-slanted) bars = 0.9o\ Based on this plot a bright- 
ness cutoff of 0.8er was selected and a minimum of seven spa- 
tial pixels was required in a valid detection. The artificially low 
(T-limit is possible because the d trous reconstruction effectively 
suppresses most of the noise. 



(10) contains the angular radius, in arcminutes, of a circu- 
lar source subtending an equivalent solid angle on the sky 
and column (11) the brightness- weighted radius measured 
from the integrated intensity map. Column (12) contains 
the solid angle subtended by the cloud in square arcmin- 
utes; (13) the number of spatial pixels in the cloud; (14) the 
number of emitting voxels in the cloud, and (15) presents 
the total integrated intensity J T m bcfo. The peak brightness 
temperature is recorded in (16) and the local RMS noise 
temperature in (17). Column (18) contains flags noting if 
the cloud touches a survey boundary in b (=X), I (=Y) or 
v (=Z). The cloud is flagged with an 'M' if multiple veloc- 
ity components were detected, i.e., if the velocity range is 
greater than the expected velocity width of a single NH3 
spectrum. Here we assume a FWHM for each line-group of 
3kms _1 , which is typical of gas forming massive stars. The 
expected velocity range of an NH3 (1,1) spectrum from an 
isolated cloud is then 44.6 km s -1 . Clouds which exhibit sig- 
nificantly different peak and brightness-weighted velocities 
(> 5kms _1 ) are likely to contain multiple sub-clouds over- 
lapping in position and are flagged with an 'E'. A small 
number of clouds detected in the NH3 (1,1) bandpass are in 
fact broad line-wings of NH3 (2,2) encroaching on the (1,1) 
bandpass and are flagged as artifacts with an 'A' in column 
18. Similar artifacts in the (2,2) data are also marked with 
an 'A' in column 18. 



5.3 Catalogue completeness 

5.3.1 Spurious sources 

Although the data have been corrected for large scale in- 
homogeneities in the noise properties by fitting polynomial 
baselines (see Section 12.21 and Appendix [Al. higher order 
fluctuations are still present and masquerade as real emission 
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Figure 12. This figure shows the percentage completeness as a 
function of peak brightness temperature. The four curves illus- 
trate the effect of varying the minimum number of spatial pixels 
required in a detection. 



near the sensitivity limits. To constrain the number of spuri- 
ous sources as a function of cloud size and sensitivity we ran 
the source finder on an 'empty' test cube at a range of sigma 
limits. The test region was drawn from the NH3 (2,2) dataset 
in the outer Galaxy (299.96° > I >292.54°, |6| <0.456°) 
and contains no believable molecular emission: i.e., when the 
source-finder is run at a very low level on both the NH3 (1,1) 
and (2,2) data, none of the detected emission is common to 
the two datasets. We would expect the J,K = (1,1) and (2,2) 
transitions to arise in the same gas and the NH3 (2,2) line 
to be weaker, hence, the detected (2,2) emission without 
corresponding (1,1) emission is likely spurious. 

Figure [TT] shows the number of spurious sources de- 
tected as a function of pixel count and sensitivity limit. 
Decreasing both the sensitivity and pixel limits has the ef- 
fect of increasing the number of spurious sources found. In 
particular, dropping the sigma-limit from 0.8a to 0.7a re- 
sults in an sharp increase from eight to seventy spurious 
detections. The 0.7a" results also show that decreasing the 
minimum size of the clouds significantly increases the num- 
ber of spurious detections. Based on these results we chose 
a sensitivity limit of 0.8<r when source finding and required 
a minimum of seven spatial pixels in a valid detection. 

Data in the test cube is of particularly high quality 
compared to the remainder of the survey and may not be 
wholly representative of all HOPS data (see Figure[7]). Fifty- 
eight spurious NH3 (2, 2) sources (i.e., without associated 
NH3 (1, 1)) were detected in the full survey, implying that 
significant numbers of spurious sources exist in the cata- 
logues. However, most of these detections have signal-to- 
noise ratios below five and subtend less than twelve pixels 
on the sky. More conservative constraints may be imposed 
when defining a high-reliability catalogue, e.g., a cutoff of 
13 pixels per source would provide a reasonable 'safety mar- 
gin' for data where the noise properties are worse than in 
the test region. Applying such a cutoff to the NH3 (1, 1) and 
(2,2) catalogues yields source counts of 523 and 198, respec- 
tively - a reduction of ~ 20 percent. 
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Table 1. Properties of the objects used to illustrate the survey 
completeness. The five objects are taken to be representative of 
an isolated low-mass core and low, intermediate and two high- 
mass star formation regions. A constant NH3/H2 abundance of 
10 was assumed. , TH and t were calculated using RADEX, 
as described in Section 15.3.31 
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Figure 13. Plot of expected main-line brightness temperature 
versus distance for the five objects in Table[T] The one, three and 
five sigma sensitivity limits are plotted as horizontal dashed lines. 
For the parameters outlined in Table [l] we should detect a 10 Mg 
cloud out to a distance of ~ 1.0 kpc, and a 400 Mq cloud out to 
3.2 kpc at the 5cr level. The turnover occurs as a result of beam 
dilution: once the source is far enough away it no longer fills the 
main beam. 

5.3.2 Completeness curves 

We have estimated the completeness limits of the catalogue 
by injecting artificial point sources into the NH3 (1,1) data 
and attempting to recover them using our duchamp pro- 
cedure. The region of Galactic plane from / = 292.4° to 
1 = 297.7° is devoid of significant emission and was used 
as an input. During the experiment the FWHM Gaus- 
sian linewidth of the injected sources was set uniformly 
to 3 kms -1 , typical of high-mass star-forming clouds. The 
distribution of peak brightness temperatures ranged from 
0.2 K to 1.3 K, spanning the expected sensitivity limits. One- 
hundred 3D Gaussian sources were inserted into the cube at 
random, non-overlapping positions, and the emission finding 
procedure was called to find them. After thirty iterations the 
positions of the clouds found by duchamp were matched 
with the injected positions and the recovered sources di- 
vided into brightness bins. Figure [12] plots the percentage 
of sources recovered as a function of peak brightness. Apart 
from the obvious sensitivity limit, the limiting parameter 
on the number of clouds found is the number of pixels al- 
lowed within a valid detection. The figure shows curves for 
Npi xe i a ^ 4, 7, 10 and 13 pixels with the la, 3a and 5a sen- 
sitivity limits overplotted. When creating the HOPS NH3 
cloud catalogue we imposed a limit of ^ 7 pixels per valid 
detection, resulting in a 60 percent completeness level at 3a 
(Tmb = 0.6 K). The catalogue is 100 percent complete at the 
5a level (T mb = 1.0 K). 



5.3.3 Completeness illustration 

To illustrate the completeness limits, we calculated the ex- 
pected NH3 (1,1) brightness temperature as a function of 
distance for five objects taken to be representative of an 
isolated low-mass core and low, intermediate and two high- 
mass star-forming regions. Starting with spherical clouds 



of a given mass and uniform density, we calculated the 
expected physical radius, projected angular size and NH3 
column density (from the H2 column density assuming 
an H2/NH3 abundance). The radiation temperature Tr, 
excitation temperature T cx and optical depth r of the 
NH3 (1,1) line was the n calculated using the r ade>Q radia- 
tive transfer package (jvan der Tak et aLll2007h . RADEX IS a 
one-dimensional non-LTE code that takes as input the ki- 
netic temperature Tkm, H2 density, NH3 column density and 
NH3 linewidth. For each molecular transition it returns val- 
ues for Tr, T ex and r. The measured main-beam brightness 
temperature T m b was calculated from Tr by applying the 
beam filling factor as a function of source distance. 

Figure [13] shows plots of expected NH3 (1,1) main- line 
brightness temperature versus distance for five objects with 
parameters outlined in Table [T] Taken at face value we 
should be able to detect sources of 1 Mq , 10 Mq , 100 Mq , 
400 M and 1000 M Q out to distances of 0.4 kpc, 1.0 kpc, 
2.2 kpc 3.2 kpc and 4.2 kpc, respectively, assuming a 5a limit. 
Dropping to a 3a level allows the detection of a 1 Mq at 
0.5 kpc and boosts the other distances by a factor of 1.3. 
This will not hold in regions like the Galactic centre where 
the bright extended emission and large linewidths effectively 
increase the detection threshold by orders of magnitude. 



5.3.4 Comparison with other surveys 

The only comparable single-dish NH3 (1,1) survey cond ucted 
in the southern hemisphere is that of lHill et al.l (2010) who 
used the Parkes radio-telescope to survey 224 high-mass 
starforming regions. The ~ 1 arcmin beam FWHM used 
means their data are four times less beam-diluted than 
HOPS and with a median ar nlb = 0.05 K the data are 4.5 
times more sensitive. As a sanity check we cross-matched the 
HOPS NH 3 (1,1) catalogue with the 104 lHill et all (|201Ch de- 
tections within the HOPS area. In total seventeen sources 
are matched within a 2 arcmin radius, accounting for the 
brightest sources in the Hill et al. catalogue. The unmatched 
source tend to be weaker and we would not expect to de- 
tect them in HOPS, within the uncertainty imposed by the 
unknown beam-dilution factor. 
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Figure 14. Plot of the 5ct detection limit distance versus mass of 
gas above the NH3 critical density (n cr i t ~ 10 4 cm -3 ) for clouds 
with representative excitation conditions. The dotted line shows 
the distance of the Galactic centre at 8.5 kpc. 



5.3.5 Galactic coverage 

The globally-averaged volume density of molecular clouds 
is typically a few 10 2 cm -3 , and the distribution of molec- 
ular gas across the Galax y is well-traced by CO emission 
from various surveys (e.g.. [Dame et ail 1200 ll . [jackson et al.l 
2006). The fraction of gas in these clouds at densities of 
10 3 - 10 4 cm , from which NH3 emission can be detected, 
is usually small - of order a few percent. However, it is in the 
high density gas that star formation occurs. Recent studies 
suggest that a critical factor controlling the star formation 
rate within a mole cular cloud is the fraction of gas above this 
density threshold IjLada et al.ll2012l ). A census of this dense 
gas fraction is clearly very important. As a large, blind sur- 
vey tracing gas at high critical density, HOPS is potentially 
a powerful tool for this purpose. 

We have used radex to estimate how the 5a NH3 (1,1) 
detection limit translates to a completeness limit for the 
mass of gas above the critical density (~ 10 4 cm~ 3 ) within 
molecular clouds across the Galaxy. Figure [14] illustrates 
the 5a distance limit as a function of dense gas mass for 
clouds with representative temperatures and linewidths. In 
the analysis we assume a fixed NH3/H2 abundance of 10 -8 . 
We expect to detect all clouds with >5,000 Mq of dense gas 
which lie between the Sun and the Galactic centre. Clouds 
with >30,000 Mq of dense gas will be detected at distances 
corresponding to the other side of the Galaxy (d w 17 kpc) . 
However, the distribution of molecular gas in the Galaxy 
is not uniform - one-third is concent rated within 3 kp c of 
the Galactic centre (see the model bv lPohl et alj [20081 ). In 
practice this means that over 90 percent of clouds which con- 
tain 20, 000 M of dense gas will be detected in the HOPS 
survey area. The HOPS survey samples approximately 66 
percent of the volume of the Galaxy where molecular clouds 
are formed. 



6 ENSEMBLE CLOUD PROPERTIES 

In the following sections we investigate the bulk properties 
of the clouds in the catalogue. 

6.1 Measured properties 

Figure [TS] presents selected plots of the measured properties 
for all NH 3 (1,1) and (2,2) clouds. 

6.1.1 Galactic longitude and latitude 

The distribution of detections as a function of Galactic 
longitude is shown in the top-left panel of Figure 1151 Be- 
tween 0° < Z < 30° the number of detections is relatively 
flat. The CMZ region, which dominated the distribution 
of integrated intensity in Figure [3] is detected as a sin- 
gle cloud, while small peaks in the distribution (e.g., at 
I ~ 13.5° and I ~ 27.5°) are due to large star- forming 
complexes visible in Figure [T] The longitude range spanning 
-30° < I < -20° contains the bright G333. 0-0.4 high- 
mass star-forming cloud and the 'Nessie' filament, giving 
rise to the peak at I w —25. There is also a small peak 
in the number of detections at / = —53°, corresponding to 
the G305 star-forming complex and a tangent point in the 
Scutum-Centarus spiral arm. The longitude distributions of 
the NH3 (1,1) and (2,2) catalogues show similar trends, with 
the NH3 (2,2) catalogue having half as many sources per bin, 
as expected from the ratio of total number of sources in each 
catalogue. There are no detections in the —60° < / < —70° 
range. 

The Galactic latitude source distribution is shown in 
the top-centre panel of Figure [15] and is significantly biased 
towards negative latitudes. At b > 0.2° the NH3 (1,1) source 
count falls from less than 30 to approximately 20 sources per 
4 arcmin bin. In contrast, no drop is seen at negative lati- 
tudes (b < —0.2°), instead the source counts remain close to 
45 per bin. As before both the NH 3 (1,1) and (2,2) distribu- 
tions exhibit similar shapes. 

6.1.2 WLSR 

The distribution of brightness-weighted i>lsr over the 
±200 km s bandpass is shown in the top-right panel of 
Figure 1151 Similar trends are seen between the NH3 (1,1) 
and (2,2) catalogues. The histogram is dominated by the 
sharp peak at wlsr ~ — 45kms _1 . Upon examination of the 
full l-v plot for the survey we see that that this feature is at- 
tributable to the sum of velocity components from the south- 
ern Galactic plane between —70° < I < —20°. Discount- 
ing the peak, the shape of the distribution is symmetric, 
with the majority of sources falling between — lOOkrns -1 < 
vlsr < 100 kms" 1 . The NH3 (1,1) source count beyond 
|vlsr| > 100 kms -1 drops off rapidly from ~ 30 to < 5. 

6.1.3 Sensitivity and angular size 

The mid-left panel of Figure [15] shows the distributions 
of signal-to-noise ratio (S/N = T mb /o"T mb ) for the NH3 (1,1) 
and (2,2) clouds. In practice the brightness and size limits 
we imposed on the source finder mean that we are sensi- 
tive to clouds with S/N^ 3, although the majority of clouds 
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Figure 15. Properties of the clouds in the NH3 (1,1) and (2,2) catalogues. The NH3 (1,1) catalogue is represented by grey histograms 
and black symbols, while the NH3 (2,2) catalogue is represented by red symbols and hatched histograms. In the lower-middle plot the 
blue hatched histogram is the subset of unresolved NH3 (1,1) clouds. See Section [6] for details. 



fall below S/Nss 20. A typical isolated cloud detected in the 
NH3 (1,1) dataset has a FWHM linewidth across the main 
group of 2kms _1 , equivalent to five spectral channels. In- 
tegrating the data across five channels would increase the 
signal-to-noise by y/5 i.e. a 3cr detection would become a 
6.7a 'integrated' detection. 

The middle panel shows the distribution of brightness- 
weighted radii r w measured from the clouds. On average 
the NH3 (2,2) clouds are smaller than the (1,1) clouds. In- 
terestingly, the NH3 (2,2) histogram exhibits a single peak, 
while the NH3 (1,1) distribution is split into two peaks of ap- 
proximately resolved and unresolved sources. We attribute 
the excess of unresolved detections to a population of spu- 
rious sources similar to the 58 low-level detections nagged 
out of the NH3 (2,2) catalogue. The thick black histogram 
illustrates the initial NH3 (2,2) detections including spurious 
sources, flagged for having no associated NH3 (1,1) emission. 
The flagged sources are clearly responsible for a second peak 
corresponding to an excess of unresolved sources. 

The mid-right panel of Figure [15] illustrates the corre- 
lation between radius and signal-to-noise ratio. Forty-four 
percent of the detections in the survey are are unresolved 



(r v < 1 arcmin) and the majority of these sources have 
S/NsJ 6. 



6.I.4 Integrated and peak intensity 

The lower-left panel of Figure fK] shows the distribution of 
integrated intensities for the two catalogues. The shapes of 
the J,K=(1,1) and (2,2) distributions are similar, although 
the NH3 (2,2) integrated intensities are weaker on average. 
Again the NH3 (1,1) distribution exhibits a twin-peaked pro- 
file, which can be accounted for by the division into resolved 
and unresolved populations. This is demonstrated in the 
lower-middle panel, which plots the J T m bdf distribution for 
unresolved NH3 (1,1) clouds overlaid on the total distribu- 
tion. The unresolved subset clearly corresponds to the less 
intense peak, which likely contains significant numbers of 
spurious sources. 

The difference in the integrated intensities is in part 
due to the higher average brightness temperature of the 
NH3 (1,1) line, as shown in the bottom-right panel of Fig- 
ure !15l The NH3 (1,1) and (2,2) distributions exhibit similar 
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Figure 16. Left panel: Distributions of galactocentric radii for the NH3 (1,1) (solid histogram) and (2,2) (hatched histogram) clouds. 
The two peaks correspond to the 5 kpc Galactic ring at positive and negative I. Middle panel: Distributions of near (left-slanted hatch) 
and far (right-slanted hatch) kinematic distances for the NH3 (1,1) clouds calculated from the brightness- weighted ■ulsr- Right panel: 
Near and far projected Galactic scale-heights for the NH3 (1,1) clouds. Note: In the first two panels negative values on the x-axis signify 
clouds in the southern Galactic plane (i.e., I < 0). Distances were not determined for the sources with \l\ < 5°. 



shapes, although the median T m b for the J,K=(1,1) transi- 
tion is 0.99 K compared to 0.87 K for (2,2). 

6.2 Kinematic distance and Galactic structure 

The Galactic disk has a well measured rotation curve which 
can be used to solve for an approximate 'kinematic' distance 
given a line-of-site velocity towards a particular Galactic 
longitude. This method has the disadvantage that sources 
located within the radius of the Sun's orbit have two valid 
distance solutions, requiring further information to distin- 
guish between the near and far distances. In addition, local 
velocity deviations due to streaming motions in Galactic spi- 
ral arms lead to a distance uncertainty of order ±1 kpc (e.g., 
iNguven Luong et al.|[20lTI ). 

Outside the CMZ the NH3 (1,1) spectra of most clouds 
exhibit uncomplicated line profiles (i.e., moderate optical 
depth with only one or two velocity components) and the 
measured vlsr is an accurate representation of the cloud 
systemic velocity. Distance estimates made using NH3 lines 
are hence more accurate than those made using optically- 
thick tracers (e.g., CO (1-0)), whose line profiles may be 
severely distorted. Bright maser lines (e.g., the 22 GHz H2O 
line) have the advantage of being detectable at further dis- 
tances, however, individual spectral components are often 
offset from the systemic veloc ity by greater than 15kms _1 
(e.g., ICaswell fc Phillipsll200Sl ). 

6.2.1 Kinematic distance 

Using the rotation curve of iBrand fc Blitzl (|l993l ) and the 
intensity weighted l-b-v coordinates we have derived the 
unique galactocentric radius, near and far kinematic dis- 
tances, and projected scaleheights for each NH3 cloud. The 
model assumes a distance to the Galactic centre of 8.5 kpc 
and a solar velocity of 220 k ms" 1 . Although recent work, 
such as by iReid et ail (2009), has suggested that the So- 
lar rotational velocity is higher than the IAU value of 
220 km s -1 , for simplicity we adopt the rotational velocity 
of the original model. A more detailed analysis incorporat- 
ing the recent developments will be presented in a future 
paper. Galactic rotation is known to depart from the ap- 
proximations of rotation curves towards the inner Galaxy, 



as such we omit sources between —5° < I < 5°. Kinematic 
distance estimates were determined for 576 clouds (86 %) 
and are recorded in the online version of the HOPS cata- 
logue. The left panel of Figure [TBI shows the distribution of 
galactocentric radii for the NH3 (1,1) and (2,2) clouds, with 
negative values indicating positions in the southern Galactic 
plane (I < 0). The NH3 (1,1) and (2,2) distributions are sim- 
ilar and we continue our analysis using the NH3 (1,1) clouds 
only, as they are more numerous. The two strong peaks at 
approximately —5.7 kpc and +4.2 kpc correspond to the so 
called 'Galactic ring' feature, a region that is most likely a 
supposition of several spiral arms, and which is thought to 
contain most of the molecular gas outside of Central M olecu- 
lar Zone l|Burton et al ] | 19751 ; IScoville fc Solomonll 19751 ). The 
difference in the galactocentric radii of the peaks may be at- 
tributed to the negatively-biased longitude coverage, but it 
is also consistent with their origin in an asymmetric struc- 
ture, such as a spiral arm. 

In the middle panel are plotted histograms of near and 
far kinematic distances for the NH3 (1,1) detections. By de- 
sign the two distributions have very little overlap, with most 
near distances falling between two and eight kpc (inside the 
tangent circle), and most far distances between eight and 
fourteen kpc. Clouds with masses or kinetic temperatures 
bright enough to be detected beyond the Galactic centre 
will be necessarily rare (see Figure I14[l , implying that the 
near distance is favoured for the majority of sources. Sup- 
port for choosing the near kinematic distances in at least 
half of the NH3 (1,1) detections comes from comparing the 
catalogue with infrared dark clouds (IRDCs). IRDCs are 
seen as dark extinction features against the bulk of dif- 
fuse Galactic infrared emission and, due to their nature, 
are t hought to lie at near kinematic distances (|Simon et al.l 
2006). We cross-matched th e NH3 (1,1) c a talogu e with the 
list of IRDCs compiled by ISimon et all {2006) from the 
8.3 um Midcourse Space Experiment (MSX) images. Fifty 
percent of the NH3 (1,1) cloud centres were found to lie 
within 2 arcmin of an IRDC extinction peak, confirming 
their association. 

We urge caution in considering the scale-height dis- 
tribution of clouds in HOPS. The b ± 0.5° coverage likely 
truncates the true distribution of molecular gas, as evi- 
denced by the flat emission profile in Figure [3] Near and 
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Figure 17. Face-on artistic rendering of the Galaxy showing the 
positions of the clouds found by HOPS at the near and far kine- 
matic distances (black crosses and gree n circles, respectively) . 
This plot is modeled after Figure 9 of lUrauhart et al.l |201ll) . 
Dotted lines mark the distance from the Galactic centre (blue) 
and the Sun (red) in kiloparsecs. The four main spiral-arms are 
labeled in black text. The position of the rotational tangent point 
is delineated by a white circle. 

far scale-heights are plotted in the right-panel of Figure 1161 
and have standard-deviations of 13.1 pc and 30.9 pc, respec- 
tively. Both values are within the expected scaleheights of 
massive star-formation tracers found in the literature, e.g., 
lUrquhart et al.l (|201ll ) measured an average scaleheight of 
29 ±0.5 pc for massive young stellar objects in the Red MSX 
source survey. 

6.2.2 Galactic structure 

Figure [TT1 illustrates the positions of the NH3 (1,1) clouds on 
an artists rendering of the Galaxy as viewed face-on. The 
imag^f] has been constructed by Robert Hurt of the Spitzer 
Science Centre in collaboration with Robert Benjamin of 
the University of Wisconsin- Whitewater and attempts to in- 
clude the most recent information on Galctic structure. The 
key features are two 'major' spiral-arms (Scutum-Centarus 
and Perseus), two minor spiral arms (Norma and Sagit- 
tarius) and nested Galactic bars: a lon g (thin) bar at an 
angle of ~ 45° to the sun-centre line Jrlammerslev et al.l 
|2000| ; iBeniamin et al.ll2005l ; ICabrera-Lavers et alj|2008l ) and 
a short (boxy /bulge) bar at an orientati on of 20 — 30° 
jBlitz fc SpergeJll99ll ; IWeiland et alj|l994h . Both near and 
far kinematic distances have been plotted as black crosses 
and green circles, respectively. For clouds within —30° < 
I < 30° the near and far distances are widely separated, 
with the majority of the far distances lying well beyond our 
3.2 kpc sensitivity limit for 400 Mq objects. If we assume 

8 http://www.spitzer.caltech.edu/images 



near distances are correct, then the cluster of sources be- 
tween 23° < I < 30° and at a distance of 6 kpc line up with 
the interse ction of the long-bar a nd the Scutum-Centarus 
spiral arm. lUrquhart et al. |201ll ') also found a large num- 
ber of high-mass stars concentrated i n this region, while 
the B UFCRAO Galactic Ring Survey |r oman-Duval et al.l 
2009) reported a large amount of molecular material in ap- 
proximately the same area. 

Clouds between —55° < I < —40° (between white 
dashed lines in Figure I17|) correlate well with the Scutum- 
Centarus spiral arm when placed at the far distance, and 
would otherwise fall between arms. 

A more thorough investigation of the near/far distance 
ambiguity will be the subject of a future paper at which time 
we will be in a position to comment on Galactic structure 
with greater authority. 



7 SUMMARY AND FUTURE WORK 

We have mapped 100 square degrees of the Galactic plane in 
the J,K = (1,1) and (2,2) transitions of NH3 at a resolution of 
2 arcmin using the Mopra radio-telescope. The survey covers 
the region —70° < / < 30°, |6| < 0.5° with a velocity range of 
±200 km s" 1 . The median sensitivity is cr Tmb = 0.20±0.06 K 
in each 0.42 km s -1 spectral channel. We have developed an 
automatic emission-finding routine based on duchamp and 
used it to find clouds of contiguous emission in the l-b-v data 
cubes. The following are our main findings: 

(i) We have detected NH3 emission across the Galactic 
longitude range between —60° < I < 30°, however, the 
outer Galaxy between —70° < I < —60° is devoid of emis- 
sion above our sensitivity limit. The Central Molecular Zone 
(CMZ, \l\ < 2°) contains 80.6 percent of the NH 3 (1,1) in- 
tegrated intensity, which appears as a single giant cloud at 
2 arcmin resolution. Most of the remaining gas is concen- 
trated within 2° < \l\ < 30°. Within the CMZ the NH 3 (1,1) 
integrated intensity peaks near the Galactic midplane, but 
the distribution with b is largely flat outside the CMZ. 

(ii) Using our emission-finding procedure we detected 669 
NH 3 (1,1) clouds and 248 NH 3 (2,2) clouds, and measured 
their basic properties. Forty-four percent of the NH 3 (1,1) 
clouds are unresolved in the 2 arcmin Mopra beam. The full 
HOPS NH 3 (1,1) catalogue likely contains significant num- 
bers of unresolved spurious sources below 5cr. However, a 
high- reliability catalogue may be constructed by restricting 
the number of pixels required in a detection to 13, reducing 
the source-count by 20 percent. 

(iii) The NH 3 catalogue contains clouds detected down to 
a 3cr level (0.6 K) where it is 60 percent complete. The cat- 
alogue is ~ 100 percent complete at the 5a level (1.0 K). As 
an illustration of the sensitivity, we would detect a typical 
400 M NH 3 (1,1) cloud (T kin = 20 K, n H2 = 10 4 cm~ 3 ) out 
to a distance of 3.2 kpc at the 5<r level. Similar clouds con- 
taining 5,000 M and 30,000 M© would be detected at the 
Galactic centre and on the far side of the Galaxy, respec- 
tively. 

(iv) We have estimated the near and far kinematic dis- 
tances towards the NH 3 (1,1) clouds. Given our sensitivity 
limits the majority of clouds likely lie at the near distance. 
Sources between —55° < I < —40° may lie at either the near 
or far distance in the Scutum-Centarus spiral arm. Further 
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investigation is required to distingush between the two pos- 
sibilities. 

Combined with other Galactic plane surveys the HOPS 
catalogues provide an invaluable tool for the investigation 
of Galactic structure and evolution. The next paper in the 
series (Longmore et al. in prep.) will present an analysis of 
the NH3 spectra, the physical and evolutionary status of the 
detected clouds, and the procedures developed to automate 
the line fitting. Further papers will publish the other spectral 
lines observed, including the J,K = (3,3) (6,6) and (9,9) NH3 
transitions, HC 3 N (3,2) and H69a. 

7.1 Data release 

The NH3 data and cloud catalogues are available to the com- 
munity from the HOPS website, http: / /www.hops.org.au[ 
via an automated cutout server and catalogue files. In addi- 
tion to the FITS format data-cubes we have made available 
the emission-finder masks, integrated intensity and peak 
temperature maps. Interested parties may also download the 
emission-finding and baseline-fitting procedures, which have 
been written in the python language. 
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Table 2. Properties of the NH 3 (1,1) clouds identified in HOPS. 
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39.2 


12.3 


60.5 


9.4 


8.1 


275.50 


1102 


18506 


2751.12 ± 10.16 


2.17 


0.17 


X 


G002.507 


- 0.032 - 


02.0 


2.505 


-0.029 


-2.4 


2.508 


-0.026 


-9.0 


-27.7 


20.9 


4.0 


1.9 


49.25 


197 


8250 


962.50 ± 7.43 


1.31 


0.19 


E 


G000.860 


- 0.067 + 


045.7 


0.807 


-0.064 


47.5 


359.875 


-0.084 


11.9 


-138.5 


146.6 


34.5 


28.7 


3748.00 


14992 


2114971 


443317.78 ± 113.19 


6.40 


0.18 


X ME 


G359.084 


+ 0.108 + 


139.3 


359.078 


0.109 


138.9 


359.042 


0.141 


130.8 


108.7 


165.4 


5.4 


3.1 


91.00 


364 


17181 


2612.32 +_ 12.95 


1.56 


0.23 


ME 


G343.423 


- 0.352 - 


27.3 


343.423 


-0.351 


-27.2 


343.425 


-0.343 


-27.3 


-46.5 


-7.3 


3.7 


2.0 


42.00 


168 


920 


182.69 ± 2.99 


1.54 


0.23 




G302.008 


+ 0.058 + 


198.1 


302.008 


0.058 


198.1 


302.008 


0.049 


197.8 


197.8 


198.2 


0.9 


0.5 


2.50 


10 


13 


1.10 ± 0.32 


0.60 


0.21 


Z 



s 



* The cloud is named for the centroid position in l-b-v space. 

T The flags in column (18) have the following meanings: X = the cloud touches a Galactic longitude boundary, Y = the cloud touches a Galactic latitude boundary, Z = the cloud 
touches a velocity boundary (±200 kms -1 ). M = multiple velocity components along the line of sight, E = an extended cloud which may contain multiple sub-clouds, A = a cloud 
identified as artefacts (e.g., invading NH3 (2,2) emission or up-turned bandpass-edges). 



Table 3. Properties of the NH 3 (2,2) clouds identified in HOPS. 



Weighted Peak 



Cloud Name* 




I 


b 


"LSR 


I 


b 


^LSR 




^max 


r c 




n 


llpix 


Hvox 


/ T mb dv 


^pcak 


Trms 


Flagst 






(°) 


n i 


1ms" 1 ) 


(°) 


(°) 


(kms" 1 ) 


(kms" 1 ) 


(') 


(') 


arcmin 2 






(Kkms" 1 ) 


(K) 


(K) 




(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


G029.995 + 0.156 


- 84.1 


29.994 


0.156 


-84.1 


29.992 


0.166 


-39.2 


-84.4 


-83.6 


0.7 


0.5 


1.75 


7 


12 


0.86 ± 0.27 


0.71 


0.18 


E 


G023.332 - 0.253 


+ 096.4 


23.333 


-0.253 


96.4 


23.267 


-0.243 


60.5 


78.0 


102.3 


2.0 


1.5 


12.25 


49 


120 


16.45 ± 0.85 


0.90 


0.18 


E 


G018.859 - 0.477 


+ 065.3 


18.860 


-0.476 


65.4 


18.883 


-0.476 


65.6 


55.8 


68.2 


2.8 


2.0 


24.00 


96 


677 


136.29 ± 2.86 


1.85 


0.26 


X 


G009.125 - 0.153 


+ 043.7 


9.125 


-0.154 


43.7 


9.117 


-0.159 


43.9 


43.5 


43.9 


0.7 


0.4 


1.75 


7 


11 


1.45 ± 0.27 


0.80 


0.19 




G000.892 - 0.057 


+ 050.7 


0.819 


-0.060 


49.7 


359.875 


-0.084 


13.2 


-128.7 


152.6 


32.6 


27.5 


3339.00 


13356 


1640704 


274605.81 ± 91.23 


4.96 


0.17 


X ME 


G338.920 - 0.483 


- 36.8 


338.922 


-0.482 


-36.9 


338.942 


-0.484 


-37.5 


-37.9 


-35.8 


1.9 


1.2 


11.50 


46 


175 


16.41 ± 1.04 


1.17 


0.18 


X 


G322.258 + 0.166 


- 171.2 


322.258 


0.166 


-171.2 


322.258 


0.166 


-171.8 


-171.8 


-170.9 


0.9 


0.5 


2.50 


10 


20 


2.02 ± 0.31 


0.67 


0.16 


z 


G300.464 - 0.498 


- 198.0 


300.464 


-0.497 


-198.0 


300.450 


-0.493 


-198.2 


-198.2 


-197.8 


0.9 


0.8 


2.75 


11 


21 


2.40 ± 0.40 


1.41 


0.21 


xz 



* , t Annotation markings have the same meanings as in Table [2] 



The HOPS NH 3 (1,1) and (2,2) catalogues 21 



400 
300 
200 
100 


300 
200 
100 



_ , , , , 1 , , I , 1 , , , I 
: Livedata baselined / 

~ / 




- I i I i 1 i i i i 1 i I i i 


i , , 


- 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
~ ASAP baselined / 


1 i i i i _ 






- ■ ■ ■ ■ 1 , ■ , ■ 1 , ■ , ■ 1 , ■ , , - 



-100 +100 

Velocity (km/s) 

Figure Al. Integrated spectrum measured from an area of emis- 
sion close to the Galactic centre. The spectrum in the top panel is 
drawn from the data baselined using only a 1 st order polynomial 
in LIVEDATA. The spectrum in the bottom panel has had further 
baselines subtracted using the ASAP fitting software. 




-100 -50 
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Figure Bl. The inset image shows an integrated intensity map 
of NH3 emission in the G333 cloud. The green contour repre- 
sents the cloud boundary. DUCHAMP detects each of the NH3 (1,1) 
spectral line components individually (top and bottom panels), 
however, the sub-clouds overlap both spatially and in velocity and 
are hence catalogued as a single merged cloud. 



APPENDIX A: SPECTRAL BASELINES 

During the initial data reduction, the livedata software was 
used to fit and subtract a first order polynomial from the 
line free channels of the spectra. In the majority of cases this 
was more than sufficient to flatten the baselines, however, 
we occasionally found that two types of baseline artifacts 
remain: 

(i) Ripples and sine-waves caused by rapidly changing 
weather conditions. 

(ii) Negative bowls adjacent to very broad spectral lines 
in the vicinity of the Galactic centre. 

Type (i) artifacts are likely to contaminate the catalogue 
with spurious sources where the peaks of the ripples rise 
above the noise. To mitigate against this problem we utilised 
the ATNF ASAlfl package to fit additional polynomials 
of order 3-5 to the data. Before performing the fit, any 
bright spectral lines were masked off using a prior duchamp- 
produced mask. This approach was entirely successful in 
eliminating these baseline artifacts. 

Broad-line sources with type (ii) artifacts tend to have 
artificially clipped boundaries and suppressed flux densi- 
ties. Figure IA1I shows examples of integrated spectra mea- 
sured close to the Galactic centre. The spectrum in the top 
panel is extracted from data processed using livedata and 
GRIDZILLA only, livedata cannot identify broad spectral 
features and so fits the wings of these lines as part of the 



http: / / svn. atnf . csiro.au /trac/asap 



emission-free baseline, suppressing the real baseline below 
zero. In order to zero the baselines we iterated around the 
following procedure until no further improvement was seen: 

• Run the emission finder to produce an emission mask. 

• Expand the cloud masks by ten percent in l-b-v . 

• Mask the emission and fit a first-order polynomial to 
the line-free channels. 

The solid line in the bottom panel of Figure IA1I shows the 
same spectrum after three iterations of fitting. All negative 
bowls were successfully eliminated from the cubes. 



APPENDIX B: RUNNING THE DUCHAMP 
EMISSION FINDER 

For detailed information on the DUCHAMP algorithm and 
settings the reader is referred to the DUCHAMP web page and 
user-guide linked therein. Here we detail how the emission 
finder was used with the HOPS NH3 data. 

At the outset the reconstructed signal-to-noise cube is 
searched one spectral plane at a time down to a global 
brightness threshold. Voxels containing emission above the 
threshold are recorded and grouped with their immediate 
neighbours. The list of detections is condensed by rigorously 
merging adjacent objects within two spatial pixels (1' = 1/2 
beam FWHM) or within the velocity spread of the NH3 
spectrum. The expected velocity limits correspond to the 
frequency difference between the main and outer satellite 
groups, plus a guard band to account for the width of the 
lines. For NH 3 (1,1) this ~ 19.7 km s" 1 plus 3.0 km s -1 (=53 
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Table Bl. Critical inputs to the DUCHAMP emission finder (version 1.1.10). 



Input Name 


Value 


Notes 


fiagAtrous 


1 OS 


Perform the o, trous wavelet reconstruction. 


reconDim 


3 


Reconstruct the data in 3D mode. 


scaleMin, Max 


1-0 


Automatically determine which wavelet scales to be included in the reconstructed image. 


snrRecon 


10 


Signal-to-noise ratio required to include a wavelet coefficient at a point in the reconstructed cube. 


filterCode 


1 


Use a B3 spline filter in the reconstruction. 


searchType 


Spatial 


Search one channel-map at a time. 


threshold 


0.8 


Search the signal-to-noise cube to a level of 0.8-cr. 


flagGrowth 


False 


Do not 'grow' clumps to a lower threshold. 


threshSpatial 


2 


Merge clumps within two spatial pixels. 


threshVelocity 


53 (68) 


Merge clumps within 53 channels (22.3 kms" 1 ), the velocity difference between the central and outer 






satellite groups in the NH3 (1,1) spectrum, plus a 3kms _1 guard band. For NH3 (2,2) the threshold 






is 68 channels (29.0 km s" 1 ). 


minChannels 


2 


Require a valid detection to span at least two channels (0.85 km s" 1 ). 


minPix 


7 


Require a detection to contain at least seven spatial pixels. 



channels) and for NH 3 (2,2) ~ 26.0kms _:L plus S.Okms -1 
(=68 channels). DUCHAMP makes no attempt to separate 
contiguous areas o f emission into comp onents in the man- 
ner of clumpfind ((Williams et al.|[l994 ) or fellwalkeeFI. 
Adjacent regions of emission in which the outer satellite lines 
overlap are merged into a single object in the resultant cat- 
alogue. Figure IBTI shows an example of an emitting region 
containing several distinct velocity components which are 
counted as a single cloud by DUCHAMP. 

When the merging is completed further criteria are ap- 
plied which filter out emission significantly smaller than the 
beam (<7 pixels), or with unusually narrow linewidths (< 
2 channels, i.e., a minimum width of 0.84 km s -1 ). Finally a 
mask cube is created which uses unique integers to identify 
the contiguous emitting regions. 

To facilitate the measurement and analysis steps the fi- 
nal mask cube is used as a template to excise each DUCHAMP 
cloud from the original HOPS NH 3 data-cubes. Small FITS 
cubes are produced containing one cloud each in addition 
to 2- and 3-D masks indicating which pixels/voxels contain 
emission. 



Ofwhm = 2 arcmin and pixel dimensions of A P i X = 0.5 ar- 
cmin, there are 12.6 pixels within the solid angle subtended 
by the beam FWHM. For weak and unresolved emission we 
may only detect the peak of the beam function, hence we set 
a threshold of seven pixels per beam. We justify our input 
criteria fully in Section r5.3.1l 



Bl Emission finder inputs 

The critical inputs for the second pass of duchamp are pre- 
sented in Table [BTI These are correct for duchamp version 
1.1.10. Before running the finder we performed the a trous 
reconstruction on the data as it confers a considerable ad- 
vantage in detecting compact and weak emission. The ideal 
reconstruction parameters are dependant on the noise prop- 
erties of the cubes and we experimented with different val- 
ues of snrRecon, scaleMax and filterCode (see Table [Blj) un- 
til the residuals contained only noise. The search thresh- 
old was based on completeness tests performed using arti- 
ficial sources injected into the data (see Section l5.3.ip . A 
0.8cr threshold was used, equivalent to ~0.16K and sensi- 
tive enough to sample the completeness curve below the one 
percent level. The last critical parameter is minimum num- 
ber of pixels required in a detection. With a beam-size of 



10 Part of the STARLINK software suite, which is available at 
http: / / starlink. j ach. hawaii . edu / starlink 



